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In  this  v/OTk  a  numerical  simulation  of  the1  mold  cavity  filling  :prdcess  was  attempted. 
Th“  midi riff! led  in  th.is  simulation  is  a  disk  which  Hot  polymer  melt  enters  through  a 
tubular  entrance  located  at  the  center  of  the  top  plate.  The  tube  is  2.54  cm.  long  and 
has  a  radius  of  0.24  cm.  The  plate  separation  and  outer  radius  of  the  disk  cavity- may 
•s  varied.  A  constant  pressure  applied1  at  the  entrance  of  the  tube  causes  the  flow. 

■:he  cavity  wails  are  kept  at  various  low  temperatures.  The  reported  results  are  for  s 
rigid  po.ly.vi nyi  chloride  (PVC)1, 

The  general  transport  equations.,  i.e.  continuity-;  momentum,  and  energy,  for  a  con¬ 
stant  density  power  jaw  fluid  are  used;  to  solve  the>  flow  problem.  It  is  assumed  that 
the  outer  radius  of  the  disk  is  very  much-  greater  than  the  plate  separation,  that  there 
is  axisymmetry,  that  only  one  of  the  viscous  terms-  in  the  momentum  equation  is  signifi¬ 
cant,  and  that  in  the  flow  direction  heat  conduction  is  negligible  compared  with  con¬ 
vention.  Constant  values  for  the  thermal  conductivity  and  Heat  capacity  of  th<?  melt 
are  used. 

Using  the  results  obtained  it  Is  possible  to  predict  "fill  times11.  The  formation  of 
q-  frozen  poly  mar  skirt  as  the  cavity  fills  may  be  followed  via  the  velocity  profiles. 

Thq  temperature  profiles  which  reflect  cooling  and  the  amount  of  viscous  heat  generated 
provide  the  basis  for  studying  resin  thermal  degradation  effects,  The  pressure  and  term 
peeature  profiles  at  the  Instant  of  fill  provide  the  initial  conditions  for  solving  the 
problems  of  packing  the  mold  and  cooling.  It  under  pressure,  from  which  final  frozen 
stains  end  residual  pressures  can,  in  principle,  be  obtained. 
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Filling  cold  mold  cavities  with  hot  polymer  melts  at  high  pressures  is  a 
complicated  physical  process.,  which  i  s  of  great  practical  interest  to  the 
chemical  and1  polymer  processing  Industries.  The  transport  approach  to  this 
process  of  solving  the  general  equations  of  change  with  suitable  equations  of 
state  to  describe  the  flowing  material  has  been  largely  ignored.  ,  No  analytic 
solution  is  possible,  and  the  non-steady  state  flow  adds  a  dimension  which 
makes  digital  computation  discouraging  because  of  the  cbre  storage  and  execu¬ 
tion  time  requirements  as  well  as  the  sophisticated  programming  required, 
these  considerations  have  forced  the  development  of  mold  filling  as  an  empiri¬ 
cal  apt  which  reties  heavily  on  experience. 

!(n.  this  work  a  numerical  simulation  of  the  mold  cavity  filling  process  was 
attempted,  The  mold;  filled  in  this  simulation  is  a  disk  which  hot  polymer 
melt  enters  through  a  tubular  entrance  located  at  the  center  of  the  top  plate. 
The  tube  is  2.54  cm,  long  arid;  has  a  radius  of  0.24  cm.  The  plate  separation 
and  outer  radius  of  the  disk  cavity  may  be  varied,  A  constant  pressure  ap¬ 
plied  at  the  entrance  of  the  tube  causes  the  flow.  The  cavity  wall's  are  kept 
at  various  iow  temperatures.  The  reported  results  are  for  rigid  polyvinyl 
chloride  (PVC) .. 

The  general  transport  equations,  Le.  continuity,  momentum,  ard  energy,, 
for  a  constant  density  power  law  fluid  are  used  to  solve  the  flow  problem.  It 
Is- assumed  tnat  the  outer  radius  of  the  disk  is  very  much  greater  than  the 
plate  separation,  that  there  is  axi symmetry,  that  only  one  of  the  viscous 
terms  in  the  momentum  equation  Is  Significant,  and  that  in  the  flow  direction 
heat  conductioh  is  negligible  compared  with  convection.  Constant  values  for 
the  thermal  conductivity  and  heat  capacity  of  the  melt  are  used. 
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The  resulting  differential  equations  are  transformed  Into  difference 
equations  explicitly,  except  for  the  energy  equation.  In  this  case  a  Six 
Point  Crank -Nicholson  Implicit  differencing  technique  was  necessary  to  assure 
thermal  stability  of  the  solution.  The  difference  equations  were  solved  by 
using  a  Fourth  Order  Runge-Kutta  integration  formula  for  the  velocity  profiles 
and  the  Thomas  method  for  the  temperature  profiles.  Convergence  to  the  dif¬ 
ferential  solutions  Is  guaranteed  but  since  a  lower  limit  was  imposed  on  the 
time  increment  by  the  core  storage  limit  of  the  computer  facilities  (27K)  and 
long  execution  times,  all  results  are  seml-quantltatlve  for  the  problem  as 
s  ta  ted . 

Using  the  results  obtained  It  Is  possible  to  predict  "fill  times".  The 
formation  of  a  frozen  polymer  skin  as  the  cavity  fills  may  be  followed  via 
the  velocity  profiles.  The  temperature  profiles  which  reflect  cooling  and 
the  amount  of  viscous  heat  generated  provide  the  basis  for  studying  resin 
thermal  degradation  effects.  The  pressure  and  temperature  profiles  at  the  In¬ 
star. t  of  fill  provide  the  Initial  conditions  for  solving  the  problems  of  pack¬ 
ing  the  fin.  1  and  cooling  it  under  pressure,  from  which  final  f rozen' sta Ins  and 
residual  pressures  can,  In  principle,  be  obtained. 

Finally,  because  so  much  of  the  total  pressure  drop  Is  dissipated  In  the 
entrance  tube,  and  so  much  viscous  heat  Is  generated  there,  this  study  In¬ 
dicates  that  the  design  of  the  gate  and  runner  system  Is  perhaps  the  most  Im¬ 
portant  facet  of  success  In  mold  filling. 
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INTRODUCTION 

For  many  years  studies  of  injection  or  transfer  molding  have  dealt  With 

the  problem  of  the  mold;  cavity  filling  process.  Most  of  these  studies  have 

been  empirical  and  experimental ,  where  "moldabil i ty"  of  a  given  resin  is 

usually  related  to  how  far  {he  resin  penetrates  a  mold>of  a  certain  shape, 

1  2 

such  as  the  spiral  mold.  *  Another  useful  empirical  criterion  used  in  prac¬ 
tice  to  insure  successful-  molding  is  that  of  having  a  resin  viscosity  in  the 
3  3  3  -1= 

range  of  10  to  2  x  10  poises  at  10  sec  at  the  molding  temperature.  Of 

course,  the  polymer  resin  has  to  be  thermally  stable  at  that  .temperature. 

As  the  industry  turned  to  molding  parts  having  more  complex  geometries, 

larger  surface  areas  and  thinner  crqsssections,  more  reliable  criteria,  had 

to  be  developed.  Thus,  recently  experimental  molds,  instrumented  with  pres- 

3  A 

sure  gauges  and  thermocouples  at  various  positions  have  been  constructed.  ’ 

No  experimental  results,  which  adequately  measure  the  pressures  and  tempera¬ 
tures  developed  |n  the  mold  as  it  is  filling,  have  been  reported  as  yet.  Fur¬ 
thermore^  and  perhaps  this  is  more  significant,,  several  people  have  made  an 

attempt  to  treat  the  mold  filling  process  as  a  problem  of  combined  transport 

3  5 

of  momentum  and  energy.  That  is,  a  "transport  phenomena"  approach  has  been 
applied  to  this  transient  and  nonisothermal  flow  problem.  The  resulting  sys¬ 
tem  of  coupled  partial  different!  a-l  equations  cannot  be  solved  analytically. 

By  making  liberal  use  of  certain  physical  assumptions  to  simplify  the  mathe¬ 
matical  problem,  a  difference  approximation  to  the  solution  may  be  obtained 
us i no  numerical  techniques  and  a  digital  computer. 

Pearson  In  his  book  "Mechanical  Principles  of  Polymer  Melt  Processing" 
has  outlined  a  numerical  method  for  solving  the  problem. **  The  mold  is  a  geo¬ 
metrically  r!mpie  flat  circular  cavity  of  uniform  shallow  depth  h.  with  an 


V 


entry  point  at  the  center.  In  deriving  the  flow  and  energy  equations  he 
assumes 

!..  Constant  resin  density  p 

2<.  Velocity  vector  is  ,given  by  ;jji(r,,x, o,oJ 

3.  Lubrication  approximation;.  <<  |~ 

.  4.  ;Axi symmetry  (no  0  dependence.^  and  symmetry,  about 
the  plane  x  «  0 

5*  Compared  with  the  viscous  term,  the  acceleration* 
inertia  and  gravitational  terms  in  the  momentum 
equation  are  negligible 

6.  Heat  conduction  Is  negligible  in  the  flow  direction 

compared  with  convection. 

7.  Viscous  energy  dissipation  (VED)  term,,  t  — •  !,s 

3u  rr  dr 

negligible  compared  with  r  0  -k—  . 

rx’  dx 

Furthermore,  he  imposes  the  following  boundary  conditions: 

a.  no  si  ip  at  the  wal  1' 

b.  there  is  a  steady  supply  of  melt  at  constant  temperar 

ture  T|  enetering  the  cavity 

c. ,_  the  beat  transfer  coefficient  L  at  the  walls,  which 

are  held  at  Tq  ,  is  constant  arid  defined  by 


'x  =  ±  h/2 


-  "  T(rV±  h/2, 


o.  the  melt  occupies  the  region  r  L  R  so  that  the  moving 
front  R(t)  is  given  by 


2  2 
R  =  a  + 


f\{i)  dt 


>''V«dr 


Using -thi$  equation  to  determine  front  position  R (t)  ignores  the  fact 
that  the  front  actual  ly  has  ,a  velocity  profile..  However,  this  I5  not  a 
serious  drawback  so  long  as  0  (t)  I's  accur~/tely  described. 

„e.  The  average  temperature  of  the  newly  formed  froa£  is 
given  by 


T [ft (t  +  At)  ,x,  t  +  Atj  = 


dx 

•  '  Q/2nR(t.\ 


The  above  equation  has  no  physical  slgnir'icance.  Averaging  over  quantities 
based  on  averages  of  known  quant  it  ies  does  not  provide  accurate  information  to. 
predict  unknowns.  Certainly  x  should  not  appear  on  the  left-hand  side. 

The  injection  system  is  further  described  by  Pearson  as  being  character¬ 
ized  :by  a  uniform  flow  rate  In  the.  early  part  of  the.  cycle  followed  by  a  uhi.r 

f  R 

form  pressure  drop,- J  dr  =  -P  ,  once  -the  pressure  drop-  reaches  P  .  This 
description  is  unsatisfactory  when  the  effects  of  the  runners  and  gate  are 
considered. 

He  then  suggests  a  numerical' method  for  solving  the  problem  as  he  deu 
scrib.cd-  It  for  a  power  law  fluid.  We  attempted  to  program  his  model,  but 
found  .certain  Inconsistencies  in  the  description  and  numerical  edition  scheme. 

it  Is  not  surprising  from  previous  comments  that  no  meaningful  average 
temperature  for  the  newly  formed  front  was  obtained  using  equation  3.  The  in¬ 
tegral  on  the  right-hand  side  of  this  equation  is  totally  insensitive  to  the 
relative  position  of  the  newly  formed  front  to  that  of  the  previous  front. 
Since  all  the  function  values  are  representative  of  an  existing  situation  and 
At  does  not  appear  either  implicitly  or  explicitly,  the  predicted  average 
temperature  is  always  Tj  .  Several  variations  of  equation  3  which  Introduced 
time  implicitly  also  failed.  We  calculated  the  detailed  temperature  profile 


---**•* 
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of  the  front  using  a  heat  balance  ;ln  that  region,  and  used  the  cud  average 
temperature  to  calculate  the  average  temperature.  We  tried  the  explicit  dif¬ 
ferencing  technique  that  Pearson  suggested;  Be  used  for  the  energy  equation 
and  .found  that  it  becomes  unstable  because  of  errors  associated  wi th  this 
approximation.  After  a  certain  point  the  predicted  temperatures  were  either 
unreasonably,  >h  I  gH  ot|  even  negati ve. 

Barrie^  describes  a  method  for  estimating  the  pressure  developed  at  the 
entry  point  of  a  similar  disk  mold  as  a*  function  of  flow  fute.  For  a  power 

law  fluid- under  isothermal  conditions,  the  integrated  momentum  equation  es¬ 
tablishes  that  the  pressure  necessary  to  support  radial  disk  flow  has  two 
components.^  One  arises  ..from  the  shear  stresses 


and  the  second  is  due  to  the  hoop  stresses 


O  #  .  r  -  1 

Hoop  stresses,  ,  cannot  be  ignored  at  high  shear  rates  (.  l(r  Ss  y  <  103sec  ), 
ou  *  * 

However,  in  the  non  isothermal  situation  a  layer  of  frozen  polymer  appears, 
in  time  which  reduces  the  effective  plate  separation.  Its  thickness  ax  cans 


be  approximately  estimated  assuming  ah  infinite  slab  of  melt  at  Tj  is  in 

contact  wi  th  a  thermaf.-ly  conducting  wa\H  at  -T  ,  provided  that  the  "freeze 

8 

off"  ;(d.e.,.  71,  or  T^)  temperature  and  thermal  di.ff.usivity  are  known.  Al¬ 
ternatively  he  proposes  a  purely  empirical  expression  to  estimate  Ax.  .  Sub¬ 


stituting  the  reduced  separation  into  the  expressions  for  the  pressure  com¬ 


ponents  and  adding  them  yields  the  pressure  estimate. 


<■  v_  ^  ^  '  •  -r-  S  '*  '  "  ‘  *  ’  -  --S  '  '-"r  '  ■*  •  ' 
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Since  the  detailed  flow -behavior -was  ignored,  It  is  not  surprising  that 
good  agreement  v/f.th  the  observed  pressures  is  obtained  only'  when  the  empirical 
estimate  of  Ax  ifs  used..  Our  results  predict  a  layer  of  frozen  polymer 
which  increases  in  thickness  as  one  follows  it  away  from  the  disk  entrance 
at  any  time.  At  the  instant  of  fill  the  difference  in  pressure  at  the  disk 
entrance  for  the  isothermal  and  non  isothermal  flows  was  only  about  k  1/2%. 

We  can,  therefore,;  confirm  that  such  a  slight  modifi cation  of  the  isothermal 
momentum  equations  should'  result  in  good  pressure  estimates.  Furthermore,  at 
high  flow  rates  the  time  required  to  fill  the  cavity  will  be  unaffected  by 
the  wal  l  temperature,-  making  the  isothermal ly  predicted  f M  I  times  reasonably. 
The  implicit  steady  state  assumption  used  in  Barrie's  work  is  the  most  mis¬ 
leading,  since  it  implies  that  the  entrance  pressure  has  a  unique  value  when 
the  cavity  is  being  filled  at  a^  constant  flow  rate.  This  is  true  at  any  in¬ 
stant  but  both  change  with  time.  To  say  tnat  the  cold  cavity  is  being  filled 
at  a  constant  flow*  rate  during  which  time  t.  a  pressure  at  the  disk  entrance 
is  also  constant,  is  inconsistent,  since  more  pressure  Is  required  to  keep  a 

colder  material  flowing  at  ,a  given  rate  than  a  hotter  one. 

3 

Harry  and  Parrot  are,  to  our  knowledge,  the  first  Investigators  to  pro¬ 
pose  and  execute  a  numerical  solution  to  the  problem  of  a  hot  melt  flowing 
into  a  cold  empty  cavity.  Choosing  a  co-ordinate  system  that  moves  along 
wl th  the  polymer-  front ,  they  simulate  the  filling  of  rectangular  slits  with 
polystyrene  at  1.1,660  psi.  The  average  velocity,  position,  temperature  and 
veloci ty ‘profiles,  of  the  melt  front  are  reported  at  various  times  during  the 
filling  process.  The  simplifying  assumptions  made  in  obtaining  the  coupled 
equations  of  momentum  and  energy  are: 


"  •'  '  '  •■ :  '  -  -  .  -  "  '  .  -  8. 

\i  The  pressure  gradient  in  the  flow  direction  Is  con-:; 

stant  at  any  given  j.nsfah.t  !•£.„*  ||:|  -s*  f(?l  „ 

'2*  The  effect -of  .cay.l.ty  edges  is  neglected. 

3.  The  fi.ui'd  i.Sr  1hc6mpresiible.._ 

k..  Conductivity,  density  and  hea't  capacity  are  constant. 


5.  -Mole!  surface  .temperature  -is  constant* 

6,  Flow  is  one  d imensdoRa 1— 

The  effect  of  the  runners  and  the  gate  are  neglected  In  the  simulation. 
They,  have  attempted  to  account  for  £ he  two  controls,  pressure  and  flow  rate, 
of  the  injection  molding  machine.  Ini  tial  1y -flow- may  be  roughly  equal  to  the 
maximum  output  rate  of  the  pump,*  since  the  hot  melt  does  not  offer  sufficient 
resistance  to  obtain  the  preset  vsour'ce  pressure.  When  sufficient  resistance 
is  experienced,  it  is.  the  maximum  source  pressure,  that  controls  the  Injection 
mold  fl 1 1 |ng  operation;  The  time  Increment  used  computationally  is  hot  preset 
but  determined  by  the  temperature  changes,  a  cons Ider&ble  advantage.  The 
ability  of  the  simulation  to  predict  short  chefs  and  successful  fills  was  ex¬ 
perimentally  verified.  However,  the  observed  and  predicted' short  shot  lengths 
or  fill:  times  were  different.--  This  is  a  worthwhile  piece  of  work.  Perhaps 
the  most  serious  assumption  -is  that  of  a  constant  pressure  gradient. 

Finally,  it  is  worth ^mentioning- a  work  which  concerns  itself  with  the 

-cooling  of  a  molding  part  after  the  cavity  has  been  filled..  The  unsteady- 
state  heat  conduction  equation  for  a  cylindrical  mold  filled  with  polyethylene 
and  cooling1  under  pressure  was  solved  numerically -by  Ken ig  and  Kama  I Their 
rigorous  .analysis  accounted  f Or  the  temperature  dependence  of  density,  specific 

heat  and  thermal  conductivity.  The  equation  o*  state  used  to  -cs cribs  the  melt 

benayioi  is 


(P  +  Tty  ’  (V  *  U*  «  R'T  . 


(6) 


9. 

The  predicted:  temperature  profiles  matched  those  measured  by  the  Instrumented 
-mold' within  ±  5°F  except  at  short  'times  in  the  region  hear  the  cold  waf'h 
What  mostly- Interested  Us  .was  the  statement  -that,  when  constant  properties 
"'ere  assumed,  the  numerical,  temperature,  stab! l  lty  depended  on  the  choice  of 
these  values  and  the  space  and1  time  increments  employed  In  the  numerical  solu¬ 
tion.  We  also  observed  that  she  choice  of  At  did  affect  thermal  stability 
and  we  did  use  constant  dens  *  ty.,  specific  heat:  and  thermal  conductivity. 

DeSCIVI.RTIOM  OF  THE  PROBLEM 

The  foil  owing  are  the  essential  features;  of  the1  problem  that  We  attempted 
tq  -oive  numerically:  A  disk-shaped  cavity,  shown  In  Fig.  1,  'the  walls  of 
~  wrjlrlr  are  kept  at  a  temperature  T  ,  is  being  filled  with  a  non-Newtonian 
polymer  melt  which  enters  the  cavity  at  a  constant  temperature  Tj  through  a 
tube  at  the  center  of  one  of  the' cavity  plates.  The  polymer  mel  t,  is  forced; 
into  thq  cavity  through  a  constant  pressure  (PRES)  at  the  tube  entrance.  Pres¬ 
sure  drops,  which  for  a  fixed  position  are  time  dependent,  occur  both  along  the 

tube  axis  and  the  disk  radial  position.  The  pressure  gradients  are  not  as- 
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sumed  to  be  l  inear  but  are -calculated. 

We  assume  throughout  this  work  that  the  melt  density  Is.  tot  a  function  of 
temperature  or  pressure.  This  is  a  serious  physical  assumption,  but  it  seems 
to  have  a  small  effect  on  the  velocity  and  temperature  fields  in  non  Isothermal 
polymer  melt  flows.  ^  The  constant-  density  assumption  dictates,  that  the  volu¬ 
metric  flow  rate  Is  a  function  of  time  only. 

The  following  two  equations  hold  for  the  physical  situations  considered 
in  this  work 


■•-■ft  Hi  f, 


Where:  PLEVEL  -  the  hydrostatic  pressure  level  at  the  cylindrical 

surface  r  *  a  . 

The  transport  equations  needed  to  solve  the  flow  and  heat  transfer  prob¬ 
lem  involved  in  the  cavity  tilii'ng  process  are  the  continuity,  r-component 
momentum  and  thermal  energy.  The  flow  is  laminar;  u  =  =  u(r,x,  t)  and 

v0  -  vx  =  0  .  There  are  no  gravitational  forces  in  the  radial  direction  and 
the  acceleration  term  is  negjected  in  the  r-momentum  equation.  We. also  further 
assume  that  u/r  and  3u/3r  <<  3u/3x,  This  is  a  good  assumption  for  cases 
where  the  me  ft  front  has  advanced  to  large  radial  positions,  but.  relatively 
poor  for  very,  short  times  in  the  cavity.  filling  process.  It  is  worth  mention¬ 
ing  at  this  point  that,  since  in  our  solution' method  we  make  use  of  the  lubri¬ 
cation  approximation  technique,  the  obtained  velocity  field  will  be,  as  stated 
above,  a  function  of  the  radial  position.  The  assumption  3u/3r  <<  8u/3x  is 
used  only  to  eliminate  a  stress  component  in  the  momentum  equation.  Thus,  the 

j 

continuity  and  momentum  equations  become: 

-  u/r  Ijr-(ru)  =  0  (9) 


1 


wheue:  S  *  l/n 


12, 


The  volumetric,  flow  rate,  assuming  symmetry  about  x  =  0  ,  is: 

n/2 


Q(t)  =  -4nr 


or 


3r* 


.  f-  i£| 
L  >rJ 
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*T)  •"* 

l/S 


dx 


W  f 
L  *-0 


G(t) 

WT 


s+i 


ax 
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(17) 


(18) 


Equation  18  ils.  used  to  calculate  (3p/8r). 

Integration  of  18  with  boundary-  condition  P^j^t)  *  0  »  Is  the  position 
of  the  melt  front,  leads  to 


p(r,t) 


;Q(t) 


.$+1 

i»7T  I 


'0;  -m 


1-n,  .  1-n 
rik  r 

1-n 


Evaluated  at  the  entrance  (r  «  a) 


p(a,t)  -  PLEVEL  « 


OT'.JBE 


»/ 


,>n/2 


s+i 


dK 


«n 


rf  "  -  a1'" 
“  T-'n 


(19) 


(20) 


The  above  equations  are  caused;  through-  m  <?  m(T).,  with  the  thermal 


energy  balance,  which  for  i.ncp^d'-iss-ible  fluids  is 


H 


PC 


DT 
P  Pt 


“  $  *  q)  ~  (v  :  Vv) 


(21 ) 


Where:  ^  —  Substantia!:  derivative 

Ui  ' 


q  a  -  KVT  the  thermal  energy  fiux 


u* 


for  the  flew  assumed  and  further  assuming  that  the  the«sisl  conductivity,  K  , 
ti  ^constant  and  that  heat  conduction  In  the  radial  direction*  (K/r^sp  fr.^f/Sr) 
}s  negtlgi b 1  e,  equation  21  becomes 


r  &T  ^  ,  3T^  ..  32T  u 

eV\3t  +  u  17/  “  K^2  “  T80  7  ‘  T 


8u 
rx  §x 


(22) 


Note  that  in  the  above  equation  the  hoop  .Stress  t00  Is  not  neglected,  whi’e 

it  is  neglected  in  the  momentum  equation.  The  Power  law  const? tut-ior*  equation 

.has  to  be  modified  because  of  the  inclusion  of  the  term  *  2  ~  . 

98  r- 


k(n-l)/2 


Thus  Tfl0  *  -  m 


A  :  A/2 1  ♦  —  becomes 

~  ~  /  r 


09 


,m 


$  [2<?>2  + 


and 


rx 


m 


#t<7>2  *  &r)n 


Thus  the  energy,  balance,  equation  22,  becomes 


PC. 


(•«  *  e)  •  >$  *  -If  *  (fl 


(23) 


where 


2m(T)|2  (^  +  {^y 


(n-l)/2 


(2*0 


The  following  boundary  and  initial  conditions  are  used: 


¥) . 

'  'X=  ± 


h/2 


“  L(T0  -  T(r,  ±  h/2,  t)) 


(25) 


where  L:  constant  host  transfer  coefficient 


T{r  =  a  »x,t)  »  T, 


T(r,x,t  *  0)  -  Tj 


The  temperature  of  the  newiy  created  front  Tfr.^.Xjtj  is  found  by  ad¬ 
justing  the  temperature  of  the  previous  front  predicted  by  the  energy  equation 
at  the  hew  time  *x.'t)  to  account  for  heat  generated  by  VED  as  well  as 

heat  lost  to  the  wails  and  empty  cavity  downstream  during  the  past  computational 
time  increment.  For  the  heat  transfer  to  the  empty  cavity  forced  convection  is 
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assumed  with  h  =  50  Kcal/m  nr  C.  The  heat  transfer  coefficient  L  is  esti¬ 
mated  by  dividing  the  thermal  conductivity  of  strainless  steel  by  a  typical 
cavity  wall  thickness,  about  1  inch.  The  value  of  L  used  is  of  the  order  of 
500  Kcai/m2hr°C. 


.  .  NUMERICAL  CALCULATING  SCHEHE 


A  numerical  simulation  of  the  flow  problem  described  is  achieved  by  re¬ 
placing  all  the  necessary  differential  and  integral  equations  by  their  finite 
differ’ence  forms.  Rather  than  represent  the  flowing  system  in  terms  of  con¬ 
tinuous  functions  of  space  and  time,  u(rjX,t)  ,,  T(r,x,t)  ,  p(r,t)  etc., 

one  uses  their  values  u,  .  , ;  T.  .  ,  ;  p.  .  at  a  o'iscrete  number  of  spatial 

•  , J  ,K  *  »J  »K  l,K 

points  and  times  given  by  the  indices  l,j,k  . 


r,  -  r,_,  +  4r, 


Xj  *  jAx  t^  *  kAt 


h  ?(L-t)Ax 


0<j  <  l 


Thus,  Uj  j  k  =*  Jnx  ,  kAt), 


The  petitions  r.  are  determined by  the  initial  condition  that  the 
mold  be  empty  at  t  *  0  and  the  boundary  condition  thatw.if  the  moving  front 
is  designated  ,  that  the  melt  occupies  the  region  r  <  f  jk  .  integrating 
the  continuity  equation  subject  to  these  constraints  we  obtain  for -constant 
volumetric  flow  rate  Q 


r2  ,  2  -t-  Q  -  at 

rlk  rik— 1  xh 


So  the  positions  r.  are  merely  previous  locations,  of  the  moving  polyner 
melt  front. 

Since  the  melt  Is  cooling  as  it  advances  into  the  cavity,  one  must  ate-? 
count  for  the  fact  that  the  polymer  wMl  freeze,  i.e. t  solidify,  when  it 
reaches  a  certain  temperature.  Secause  of  the  cold  cavity  walls  the  ces¬ 
sation  of  flow  will  start  there  and  gradually  move  inward  toward  the  center 
of  the  flow  region.  This  frozen  material  near  the  wal Is  prevents  further 
flow  ir.  this  region.  Consequently,  the  hot  polymer  finds  the  cavity  be¬ 
coming  narrower  as  cooling  progresses.  The  spacing  h  in  the  continuity 

equation  is  thus  replaced  by  h  .  AM  the  material  in  the  region 

c 

h  a£  3c-. <  h  is  frozen^ whereas,  if  0<x  <  h  the  material  is  flowing.  The 
e  “ 

position  of  the  new  front  Is  now  given  by 


The  equation  of  motion  for  a  power  law  fluid,  equation  is  numerically 

represented  In  the  following  manner.  The  left-hand  side,  PPDR(I.K).  5s 

evaluated  for  each  radial  position  and  time  by  using  the  finite  difference 

form  of  equation  18,  Therefore,  the  pressure  gradient  is  not  assumed  con- 
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stant.  Simpson’s  rule  is  used  to  approximate  the  irtegral. 


16. 


-  V,  <~* 


The  local  fluid  velocity  at  esch  height, (J)  position, is  obtained  by  ba- 
ginning  at  the  wall  (where  it  is  always  ?ero  because  of  the  no  si  ip.  boundary 
condition)  and  calculated  toward  the  centerline  using  &PDft{},K)  and  the 

I  K 

Fourth  Order  Runge-Kutta  rcheme  to  calculate  the  unknown  quantity  D'ELV, 
shown  schema  ri cal ly  below; 


DEl-V  -H  t+~ 


*  To  determine  whether  DPDR(I,K)  is  correct,  the  pressure  at  the  entrance 
of  the  disk;  PLEVEL,must  alSo  be  compuLed  from  DPDR(i,K).  If  PLEVEL 

s 

(equation,  20)  added  to  the  pressure  drop  through  the  tube  DPTUBE  equals  the 
total  pressure  drop  ?RE$,  then  the  equation  of  motion.  Is  satisfied.  If  not,. 
Q.TUBE  or.  the  right-hand  side  of  equation  18  is  appropriately  varied  until 
equation  8  holds. 

The  energy  equation,  23  and.  2k,  Is  a  parabolic,  non-Homogeriepus  partial 
differential  equation.  It  is  transformed  by  the  Crank-Nickolson  "six  point 
implicit  form"  method^  to  a  system  of  (L-1)  simultaneous  equations. 
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7«,j+*,kHrl 


+  T. 


+  T> 


W.fc 


T.  .  ,  . 


r 

L  ri'rM  J 

(30) 


where  £  *  'kAt/pCp(Ax)  is  a  dimensionless  group. 

Each  equation  is  rearranged  with  the  unknown  quantities  on  the  left  and:  cast 
in  the  matrix  form: 
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The  above  system  of  equations,  which  now  represents  the  energy  balance  in 
difference  form,  is  solved  by  the  method  of  thomas.^  Details  of  the  solution 
as  well  as  copies  of  the  computer  program  may  be  obtained  through  the  authors. 
Conduction  in  the  flow  direction  is 


k  a 
7  W 


(32)- 


I  n  d  i  f ference  fore: 


qd 


;  rrfcoH  _  * 
TROT-+  jUi-DJ/z  *' 


T(1,J,K)  -  T(l-l;j,K) 
'  ftlilr-  Rir-lj  - 


thcoh  *  Hi±LM)  -  zr(MvjO  - t  *Mh,j,k) 

tR(f)  ;  RO-tjj2 


(33) 


Convection  In- the- flow  direction  is 


pC  u 
P 


37 

3r 


» 


‘’c 


(34) 


In.  di f ference  form 


„c  .  RHGCP  *  ■»■%>.*>.  *  ..«,L(lrV-,K)  *  H;j:j)J.K)  (35 


The  CGS  units  of  either  term  are  [cal  cm  ^  sec 

Thus,  one  can  calculate  each  quantity,  approximately,  at  any  point; 

Consider  a  0.635  cm  disk  being  filled  under  an  applied  pressure  of 
7500  psi  for  a  melt  entering  at  202°C  and  oollng  because  of  30°C  cavity 
Walls. 


At  I  «  11,.  J  =  1!  and  t  -  2.75  sec,  In  the  hot  melt  region,  we  find  that 
qd  “  7.38  x  10  ^  and  q£  *  121.0, x  10  Thus,  radial  conduction  Is 
negligible  In  the  "hot  melt  region,- 

Wher.  the  rigid  PVC  in.  the  vicinity  of  the  cold  wall  cools  below  T  ,  a 
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new  sol  id  phase  is  formed  and  the  boundary  condition  (25)  changes.  Heat  from 
the  melt  ls  first  conducted  through  this  frozen  layer  and  then  out  through 
the  cold  steel  walls.  More,  correct  boundary  conditions  are 


at 

3X 


h  72 


K(sol 

d'i'soV 


j-j  T.(r,he/2,t)  -  T(r,x,t)]  (36) 


jx  ft  h/2 


-  T (r, h/2.  t) 


A  slight  change  in  the  results  is  obtained  and  since  Incorporation  of  .the 
above  feature-  for  the  * ‘frozen  skin"  region  requires  lengthy  programming,  addic¬ 
tions,  i-t  was  not  -included  Tn  the  numerical  solution  of  the  energy  equation. 
Finally,  temperature  dependence  of  thermal  conductivity  and  heat  capacity  can 
be  included  in  the  solution. 


RESULTS  _AfiD=  DISCUSSION 


Rigid  PVC  melts  were  considered  wi th  Power  Law  constants  of  n  »  0.5 
and  m  (202^0)  *  k  x  10?  poises.  The  activation  energy  for  flow  Was  taken  as 
27.6  Kcal/g  mole.  On  Fig.  2  pressure  of  the  non- isothermal  as. well  as  iso¬ 
thermal  flow  is  plotted  as  a  function  of  radius  at  various  times  of  flow.. 

For  the  non- isothermal  case  (Tj  ■  202°C  and  TQ  *  30°C).  where  both  VED  and 
cooling  are  present,  the  pressure  profile  along  the  disk  at  a  given  instant 
shows  considerable  deviation  from  the  Idealized  Isothermal  pressure  profile. 
Initially,  VED  predominates  and  lowers  the  melt  viscosity.  This  prevents  the 
pressure  at  the  disk  entrance  from  reaching,  as  high  a  level  as  It  did  when  no 

•tfC  n.  .TLa  -MAi.  1.*.-  t-—.  J  t.n«.  TL__ 

ttm«  j/t  oo\mvi  mo  iiuv  rcou  i  t  10  a  inure  yrautiar  ^rcosui  c^yrau  itsilU'  i  Rcn-f 

because  of  the  longer  contact  times  with  the  cold  walls,  cooling  becomes  con1? 
trolling.  This  Increases  the  viscosity  and  PLEVEL  rises  abdve  the  value  It 
had  when  cooling  was  neglected.  The  resulting  steeper  pressure  gradient  Is 
mitigated  only  when  the  presence  of  frozen  polymer  near  the  front  constricts 
the  flowing  melt  and  Intensifies  VED." 

On  Fig,  3  we  see  the  effect  of  varying  the  applied  pressure  on  the  nor¬ 
malized  flow  rato  (Q^/Q^)  anc*  front  position  (R/b)  for  an  0.3175  cm  thick 


dJ&fc  viffh  walls  at  3C°C.  The  Initial  /polymer  temperature  is,  2Q2°C. 

-  * r  At  5Q^&  p gi  $be  flow  *ate  decline  is  gradual  ib'ut  wi  th  a  steeper 
Imtiai-lvv  The  jfSQQ  p£i  curve  begins  above,  but  then  fails  below  the  5000 
psi  curve  by  not  losing  its  steeper  sjope  as  quickly,.  At  11,600  psi  the  flow 
curve  is  agaih  initially- above  but  ui timately  below  both  curves  representing 
the;  lower  pressures.  the  increase  viscous  heating  of  the  flows  at  the  higher 
pressures  is  responsible  for  their  obiTity  to  maintain  higher  flow  rates 
ini'fciaKiy,  because 'of  the  accompanying  drop  in  viscosity.  The  hotter  melt 
simply  has  a  lower  viscos ity.  However,  when  cooling  begins  to  predominate, 
flow  retardation  if  more  Intense  because- thii  non? isothermal  situation  is  more 
extreme.  Naturally,  increasing  the  applied  pressure  advances  the  front  more 
quickly. 

In  PVC  whenever  flow  is  accompanied  by  a  high  amount  of  VED,  the  problem 

a 

of  thermal  degradation  appears..  Since  increasing  the  pressure  has  the  ef¬ 
fect  of  raising  the  average  temperature  through  viscous  heat  generation,  more 
degradation  will  be  observed  at  higher  pressures.  The  effects  or  degree  cf 
thermal  degradation  have  not  been  introduced  in  the  numerical  solution. 

On  fig.  k  one  can  see  the- effect  of 'Varying  the  plate  separation  on  the 
normalized  flow  rate  for  a  disk  being  filled  with  a  PVC  melt  at  2C>2^C  under 
5000  psi.  The -Wall  temperature  is  ,30°C.  As  the  mold  cavity  becomes  thinner 
the  flow  rate  falls  off  more  rapidly.  The  narrower  gap'  Imposes  a  physical 
resistance  to  flpw  which  can  be  seen  In  terms  of  the  steeper  pressure  gradient 
and  higher  shear  rate.  More  energy  Is  absorbed  as  viscous  heat  generation. 
Therefore,  thermal  degradation  of  PVC  may  be  a  problem  :ln  thin-  cavl ties.  Re¬ 
sistance  to  flow  due  to  cooling  Is  also  playing  a  role  In  the  diminishing 
flow  rata.  * 

On  Fig.  5  ono  observes  the  effect  of  varying  plate  separation  on  the 
front  position  for  the  same  conditions  as  the  previous  figure.  As  expected 
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after  a  given  time  interval  the  thinnest  disk  has  the- 'largest  radius  Even 

~  .J 

with  a  relatively  high  flow  rate  sufficient  material  does  not  flaw  into  the 
thicker  disk  in  the  same  time  interval  to  reach  the  same  cadius. 

On  Fig.  6  we  observe  the  effect  of  varying  the  mejt  temperature  Tj  on 
the  normalized  flow  rate  and  the  front  position  for  a  disk  having.  0.3175  cm 
piate  separation  which  is  being  filled  under  a  constant  applied  pressure  of 
7500  psi.  When  the  wall  is  at  30°C,  the  flow  rate  of  the  coolest  melt  at 

124°C  falls  off  with  the  least  rate,  whereas  the  flow  rates  of  the  hotter 

melts  at  202°C  and  210°C  ultimately  show  a  more  drastic  decline  as  the  dif¬ 
ference  between  the’  initial  melt  temperature  and  :the  mold  temperature  in¬ 
creases.  We  notice  that  initially,  the  hotter  melts -maintain  a  nearly  constant 
flow  rate  for  a  length  of  time  which  increases  with  increasing  melt  tempera¬ 
ture  Tj  .  This  is  a  VED  effect;  the  heat  generated  at  the  nigh  Flow  rates 
lowers  the  local  viscosity  near  the  entrance  which  allows  the  polymer  to  flow 

for  a  short  time  at  nearly  a  constant  rate.  Then  cooling  from  the  30°C  walls 
becomes  controlling  and  a  drastic  drop  in  flow  rate  occurs.  In  the  case  of 

the  1 9^°C  melt  not  enough  heat  is  generated  to  make  this  effect  apparent  £nd 

only  a  gradual  decrease  in  flow  rate  is  evident.  However,,  even1  Ihg’lVEj): 

is  . apparent  when  the  Initial  portion  of  this  curve  is  compared  to  ap- assumed; . 

isothermal  situation  (no  VED' and  Tq  =  20 2° C  =  Tj)  <  The  greater  amounts  of 

VED  generated  by  the  hotter  melts  will  increase  the  amount  et’  th&rma.l  de* 

gradation  in  these  non- Isothermal  flows. 

As  expected  the  hotter  the  melt  the  further  the  polymer  front  advances 

in  a  specified  length  of  time.  For  example,  at  0.-5625  seconds  the  melt  which. 

entered  at  2lO°C  has  filled  57.7$  of  the  disk  cavity  as  compared  to  22.81  for 

that  which  entered’ at  J94°C. 

In  view  pF  Jfho  marked  difference  in  the  flow  rate  curves  for  the  nofr 
isotherma!  and  isothermal  curves  for  the  202^0  melt,  the  facr  that  thp  front 
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positions.  do  not  differ  markedly  also,  is  somewhat  surpf  _lsins»  Bc*th  cavi  ties 
arer36' 2*  ful  l  afiter  0-562$' S-ecdndSi  By  visual  inspection  it  appears  that  _ 
the  area  under  the  two  corresponding  $.  vs  t  curves  are  about  equaJ ;  it  is 
this  quantity- (volume)  that  determines  how. much  of  the  cavity  has  filled  with 
melt.  The  implication  here  is  that  as  long  as  estimating  fill  times  is  the 
primary  concerns  the  isothermal1  mode}  Is  adequate.  Simulations  under  a  number 
•.of  different  conditions  support  this  conclusion.  An  explanation  and  some  com- 
ments  regarding  its  reliability  can* be  found  in  the  discussion  of  Fig.  7- 
.  Barrie^  in  his  treatment  did  assume  that  the  front  position  R  was  in¬ 
dependent  of  the  thermal  environment..  To -account  for  hon- Isothermal  effects 
he  only  changed  h  in  the:  isothermal  equations  h  and;  5* 

* 

Since  our  work  endorses  this  assumption  and  predicts  that  the  difference 

in  PLEVEL  in  the  two  cases  is  about  k  ]/2%  (see  Fig.  8) , -Barrie's  success-  in 

estimating  PLEVEL  based  on  empirical ly  fitting  a  simple  theoretical  model  is 

-  • 

probably  more  than  a  happy  coincidence. 

Two  sets  of  physical  conditions  are  compared  on  Fig.  7.  An  Isothermal 
fill,  with  both  the  polymer  and  cavity  wall  at  202CC»  Is  contrasted  with  the 
same  0.635  cm  thick  disk  be  trig  filled  at  the  Same  pressure  but  with  the  wail 
at  30°C.  The  appl led  pressure  Is  varied  and  the  non-isothermal  and  isothermal 
f  low  processes  are  looked  at  agalri.  Increasing  the  appl  led  pressure  shortens 
the  time  required  to  fill  the  cavity,  Irrespective  of  the  thermal  environment. 
At  1-1,600  osi,  whether  the  wall  is  at  30°C  or  202°C  has  no  discernible  effect 
on  where  the  front  is  at  a  particular  instant  or  the  time  it  takes  to  fill 
the  cavity.  However,  at  the  lower  pressures  the  simulation  shoWs  that  the 
non- Isothermal  filling  actually  occurs  sooner,  although  initially  no  effect 
on  the  position  of  the  polymer  front  is  evident.  One  explanation  for  this 
effect  is  that  VED  Is  more  pronounced,  Thus,  increasing  the  applied  pressure 
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makes,  the  Injection  cycle  more  ir  sensitive  to  temperature  variations. 

in  figs.  3  and  9  a  comparison  of  an  isothermal  (T}  =  202°C  3  Tq)  and 
a  non- isothermal  (Tj"  =•  -202°C;  Tq  =  30?C)  cavity  f  Li  l  ing  process  is  made  in 
terms  of  the  pressure  developed  at  ;fhe  disk  entrance,  PbEVEL^  Results  are  il¬ 
lustrated  for  constant  applied  pressures  of  .5.000  *•  75OOandll,6O0psr.  The 
plate  separation  Is  <0.6|5  cm.  In  the  hpn- isothermal  sft.atiohs  PLEVEL  gradtr- 
ally  increases  as  material  flows  into  the  disk  cavity  apd  reaches  about  \M.S% 
of  the  applied  pressure  (PRES)‘  just  prior  to  fiM.  The  isothermal  cases  are 
characterized  by  PLEVEL  initially  increasing  more  rapidly  but  only  reaching  | 

Ilf.  ]%  of  PRES  just  prior  to  fill. 

-  *-  ~  _  a 

The  crossover  of  the  non- isothermal  and  isothermal  curves  was  also  ob¬ 
served  with  the  flow  rate  curves  (see  Fig.  3).  The  fact  that  they  occur  at 
almost  the  same  time,  .strongly;  suggests,  a  cause-effect  relationship  between 
PLEVEL  and  flow  rate,  but  it  could  also  be  fortuitous.  -Initially  the  VED  In 
the  non- Isothermal  f  ill  ihg  lowers  the  viscosity  and  prevents  PLEVEL  from 
'building  up.  As  cooling  becomes  contr.ol-Hng  PLEVEL  rises  sharply  followed  by 
a  gradual  increase  as  the  cooling  rate  diminishes.  The  absence  of  VED  in  the 
Isothermal  case  allows  the  melt  to  experience  cool  ing.  immediately;;.  PLEVEL  is 
proportional  to  the  increasing  viscosity.  As  the  pressure  is  Increased,  the 
higher  VED  generated  makes  this  difference  in  behavior  more  obvious. 

Illustrated  in  Fig.  10  for  comparison  are  the  velocity  and  temperature 
profiles  of  an  Isothermal  and  non- isothermal  moving  front  the.  Instant  Just  be¬ 
fore  ft  hi  occurs  in  a  0.635  cm  cavity  under  an  applied  pressure  of  7500  psi. 

The  presence  of  a  cold.  (30°C)  wall  drastically  decreases/  the  yelocity  gradi¬ 
ents  near  it.  Thus,  high  velocity  gradients  with  VED  predominating  ore  ob¬ 
served.  These  effefcts  are  responsible  for  the  "nipple-like  shape"  of  .the  non- 
isothermal  velocity  profile.  The  almost  stationary  polymer  near  the  cold  wall 
is  what  is  referred  to  as  the  "frozen  polymer  wail"  which  reduces  the  available 


cross-section thereby,  raising  the  shear  rate  ay, C  increasing  the  W?D  in  the 
moving  golyme'fe.  The;no'n.-iso, thermal  temperature  profile1  can  .also  be  described 
Ms a.  fia ttened  parabola  w ith  the  higher  temperatures  hear  the  center!  line. 

The  fact,  that  it  Tics  entirely  .be.] ow  the  IsCthe rmal  t  «mre r*<  t  i  ”  ,  profile;  is 
peculiar  tp  the.  case  at  hand.  Were  the.  Initial  malt  teir.pffptufs>  greater  of 
;fhe  gay;j.'.ty  more-narrow  or  the.-  applied*- .pressure  higher,  It  fs  likely  that'  some 
of  the-  center  -pprtior  would  Tie  above  the  horizontal  line  Which  represents 
the  Isothermal  temperature  profile. 

.'•Qr  Rg.  ir  we -see  the  shear  rate  profiles  corresponding  to  the  velocity 

profirfes  of  Fig;.  10,  For  isothermal  flow  the  shear  rate  Increases  from  zero 

at  the  center  line  to  a  maximum  at  the:  wall..  This  maximum  Ties  between  the 

certteriJne  and- the  wail  Ip  the  isothermal’  case.  The  cooling  effect  of  the 

cp’d  -.wail  renders  the,  polymer  in  its  vicinity  nearly  Imnipbl le,  assuring  iow- 

shear  rates  in  that  region.  Too  far  from  the  wall  tp  experience  cooling,  the 

gradually,  rising;  sheaf  rate  Is  augmented  by,  the  generation  of  viscous  heat 

which  lowers  the  viscosity.  The  sharp  rise  is  foil  owed  by  an  even  sharper 

drop  when  ’the  walls  are  near  enough  for  cooling  tb  ’become  controlling. 

On  Fig.  12  the  effects  of  the  cpmpiitstioWl  times  increment,  At  ,  on 

the  average  melt  temperature  and  the  position  of  the  front  are  shown  under 

conditions  that  result  in  fapld  fi ilTng.  The  front  position  Is  independent 
of  the  At  •employed,  but  the  average  melt  temperature;  depends  on  i  t  strongly 

For  reasons  of  compute f  storage  It  has  not  always  been  possible  to  keep  At 
properly  small  and  still  carry  the  calculation  far  enough  to  fill  cavities. 
Fortunately,  the  cool Ing  effect  is  so  overwhelming  as.  the  melt  advances  Into 
the  cavity  that  the  average  front  temperatures  based  on  different  At's  ap¬ 
proach  each  other  near  the  outer  radius  of  the  disk.  Nevertheless,  if  one 
wants  to.  calculate  the  extent  of  res  la  degradation,  due  to  Increased  tempera- 

mlt 

tures,  then  very  small  Increments  of  time,  about  2  x  10  sec,  would  have  to 
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be  used  in  the  numerical  solution. 

On  Fig,  1$  the  average  melt  temperature  as,  the  filling  process  proceeds 

t  . 

is  plotted  for  two,  pressures  and  two  wal  l  temperatures  at  At  -  0y05  sec.  ’ll' 
Is  observed  that  the  Wail  temperature  dons  not  have  nearly  as  large  an  effect 
on  the  temperature  as  that  of  pressure.  Thus,  i f  degradation  of  the  PVC  resin 
is  a  problem  lowering  the  pressure  and:  hot  .the  mold  temperature  would  be  of 
help. 

In  practice  the  ambient  temperature  of  the  mold  surface  varies  consider* 
ably.  The  thermal  environment  of  the  first  few  pieces  is  entirely  different 
than  that  of  subsequent  cycles,  since  the  room  temperature  varies  and  cold, 
mold  surface  heats  up.  Nonetheless,  it  is  not  usual  to  find  elaborate  cooling 
equipment  because  no  Significant  improvement  in  quality  or  cycle  time  isnir*' 
served.  This  5s  as  predicted  by  the -  simulation. 

In  Fig.  Ill  one  can  observe  the  effect  of  VED  in  the  tube  on  flow  race  and 
front  position.  Q  is  32.4  cc/sec  fpr  both  cases.  The  initial  polymer  tem¬ 
perature  is  202°C*  but  when:  VEp  in  the  tube  considered,,  the  melt  Irritlalty 
enters  the disk  at  •222.-4l0C*,,  The  flow  rate  declines  more  rapidly  than, 
wh*  n.  an  Isothermc!  tube  was  assumed  and  Js  accompanied' by  a  slower  rate,  of 
’filling.  The  primary  reasons  for  <he  slow  ill-l  ing  are,  that  flow  in  the  disk 
is  determined  b.v  PLEVEl,  which  the,  hotte--  melt  prevents  from  building  up,  and; 
that  the  applied  pressure  needed  tp' establish  Is  2000  psi  less  if  the 
polymer  heats  up  to  222.141°;  during,  its  travel  through  the  tube,  sa  that  thp.rp 
Is  less-  applied  to  begin  with.  It  is  also  true  that  the  velocity  gradients 
Imposed  Initially  oh  the  two  flows- will  be  more  severe  for  the  hotter  melt 
near  the  wall. 
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CONCLUSIONS 


f;  Thefollpwlng:  conclusions  can  be  reached 'from  the  numerical  solution  of 

&>*  \  '  .  ' 

p  K  the  complied  mci!?ehtym- and  energy  balance  equations  which  describe  fhe>  mold 

i  -  fill  ling  operation  ;o.f  niqld  PVC  mel  t  in  a  di:sk  caV'l  ty. 

ir '  *-■  ‘ ,  ’  _  .  : 

\  In  the  rapge  of  process  variables  that  we  investigated  ft  was  found,  that 

fey  >  •  -  *  ‘  - 

r  the  ’If  r!  1  times"  or  non^ isothermal  and  isothermal  jpelt  flow  were  almost -idanti 


cal,  -despite  the  fact  that  they’ differed. in  al 1  other  respects,  Although  we 
have  no  clear  explanations  for  this,  result,  jts  practical  implications  are 
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cleef.  There  is  np  computational  need  of  Including  the  energy  equation  in  the 
solution  of  the; /problem,  if  obtaining  mold  filling  times  is  the  only  Informa¬ 
tion  sought..  Any  further  work  on  this  problem  should  Investigate  a  larger 

j 

range  of  processing  variables  and  material  parameters  to  determine  whether 
ispfherjal  and  non- iso  thermal'  fill  times  are  nearly  equal  in  general.  The 
above  result  also  justifies  the;  isothermal  assumption  of  Barrie  who  in  his. 
treatment  attempts  to  predict  the  macroscopic  aspects  of  .the  mold  process. 

The  agreement  between  hi's  largely  empirical  treatment  and  experimental  results 
can  now  be  understood.  Furthermore,  because  he  has  looked  at  a  different 
polymer  resin, an  ethylene-propylene  copolymer  wc  suspect  that  our  findings  on 
the  Similarity  between  isothermal  and  non- isothermal  fill  are  more  general. 

The  similarity  of  Figs.  3  an  d  6  establishes  the  existence  of  a  pressure- 
melt  temperature  equivalence.  Raising  the  source  pressure  has  the  same  effect 
on  flew  rate  and  front  position  as  increasing  the  initial  melt  temperature; 
this,  s  1ml  larky-  can  be  extended  to  the  temperature  prof  i  lev 

The  full  non- isothermal  -non-Newtonian  me! t  description  of  the  problem  Is 
needed  to  obtain  detailed  Information  on  the  process f  such  as  information 
leading  to  the  establishment  of  the  extent  of  resin  viscous  heating  and' de¬ 
gradation,  distribution  of  the  residual  or  ’’frozen11  'strains  and  distribution 
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of  the  residual  pressure  <?1ong  the  radius  of  the  mold  cavity.  This  is  pre¬ 
cisely  the  information  that  is  needed  for  large  pa^t  molding  snd.:  iliac  our 

-  O  -  _ _ _ _ 

solutions  provides  approximately;  The  results  are  approximate  and  correct 
only  ic  a  semi-quantitative  sense  because  of  the  limitations  that  the 

-computpr -'tyg.iem  of  Stevens  Institute  imposed  on  the  solution..  Specie 
fical  ly,  theccmputaC’c^aldiffrceli  iescan  be  stated^as  follows:  The  prob¬ 
lem  i s  three- d i n^nsiona l frj  the  difference  form  of  the  transport  eq^ailcri 
a  time  increment  Afc  has- -to' be  used;  This  should  be  Kept  appropriately 
•small'  so  -tefrat  not  oT£y  Convergence  td-the  solution  of  the  differential 
equation  is.  ichTeved,v.-bpt"  fromvatnofe  practical  point  of  viewy  no  tempera¬ 
ture  fhstabi  jlt| es  are-obtained  numerically.  Whereas,.  It  was  always  possl- 

j  - 

ble  to  keep  At  small  enough  to  achieve  the  latter-,  the  system's  maximum 
core  capacity  of  27&  imposed  a  lower  limit  on  At  if  the  outer  radius  of 
the  mold  cavity  to  be  filled  satisfied  the  condition -that  b  >>  h  the  - 
basis  of  the  1  ubr i cation  approx imati oh.  The  execubon  ti^  for- Case  IV 
simulation  when  27Kwas  Util i zed  Was  about  28  rut  hutes.  Because  of  these 
-system  limitations  there  is  no  guarantee  that  convergence  to  the  differed? 
tlal  solution  was  achieved.  Therefore^  the  values  obtained  for  T(r,x,t)  , 
u(r-yX,t);  etc.  may.  hot  be  quantitatively  ^correct,  which  is  why  we  report 
the  .detailed  results  in  a  sem Irquenti tat iys  fashion. 

The  slpiulatJou  affords  s-  means,  -cjr  looking  at  both  the  average  ternperatui* 
and  the  detailed  temporatute  profile  at  a  particular  flow  radius.  Since  5 

w  S  ^  r  -  _ 

significant  amount  of  frozen  polymer  .may  be  present  even  v/bnn  the  .average 
temperature  is  wel  l  above  ,<?r  general  Izations  based  jiri  average  itsm? 

perature  may  be  misleading,  A  detailed  examination  of  T(rpf,.tj  on  the  Ot-he 

j  S 

hand  does  provide  insight  Into  several  phenomena  with  which' the  practical 
molder  must  contend.  For  example,  tbs  build  Uj  of  a  frozen  polymer  skjn  near 
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the  w»H  as  the  hot  melt  flows  into  the  cavSty,  and  the  creation  of  intense 
VEO  at  the  core  which  may  initiate  or  intensify  thermal  degradation  in  a  resin 
speh  as  PVC.  By  providing  the  temperature  and  pressure  profiles  at  the  in¬ 
stant  of  fill,  realistic  initial  conditions  ere  now  available- for  solving  the 
related  problems  of  hew  the  pressure  in  the  disk  increases  during  the  packing 
stage  and  how  the  molded  part  cools  under  pressure.  From  these -solutions 
final  frozen  strains  in:  the  material  and  the  residual  pressure  may  be  obtained 
Both  are  intimately  related  to  the  quality  of  the  molded  part. 

For  r*ori-isothefi.io:I  ^lows  which  ignored  VED  in  the  tube  entrance,  the 
.pressure  at  the  disk  entrqnde  was  \k.$%  of  the  source  pressure  at  the  instant 
of  fljt.  Under  these  condi Cions  the  gate  and  runners  have,  a  profound  effect 
on  lower i ng  the  efficiency  of  the  cycle  by  dissipating  the  bulk  of  the  applied 
pressure  before  it  has  a  chance  to  advance  the  melt  into  the  cavity.  If  r.his 
Were  available  at  the  disk  entrance  there  is  ho  question  that  the  mold  would 
ffHmiich  faster.  Melt  which  is  heated  by  VED  during  its  travel  throuor  th* 
tube  will' -be.  delivered  at  the  disk  entrance  either  at  the  same  flow  rate  under 
a  lower  source  pressure,,  or  at  a  higher  flow  rate  under  the  same  pressure, 
the  development  of  PLEVEl  will  b.e  hindered  by  the  lower  viscosity  of  the  melt 
•exiting  from  the  tube  and  it  should  not  reach  llf.b*  at  .the  fill  instant  in 
either  case-  Consequently,  5  icw  PLEVEL  and  q  lower  source  .pressure  might 
result  in  longer  -fill  cisnes  because  of  the.  increased  flow  rate;.  This  will 
general?*  more  intense  viscous  healing  throughout  the  cavity  as  it  fills. 

Based  on  these  observations  We  conclude  that  the  gate  and  runners  are  the  most 
important  factors  to  ae  considered  in  mold  design-  1'gTiorlng  their  effect  in 
future  studies  may  be  the  most  important  reason  for  falling  to  predict  actual 
fl-Tllng'  behavior. 

It  is  not  ^urfirising,  2-n  view  of  this  study,  why  the  fancy  design,  of  the 
runner  system  has  been  responsible  for  successfully  f ill  ing. large  radios,  By 


. 

carefuSly  designing  the  corners  in  runners  and  using  hot  rsultigated  runner 
manifolds,  the  adverse  pressure  find  viscous  energy  dissipation  effects  of  the 
entrance  can  be  ;iessended. 

Experiments  to  quantitatively  support  or  refute  the  predicted  flow  be¬ 
havior  would  be  the  most  valuable  supplement  to  this  work.  Furthermore,  the 
work  should  be  completed  for  rigid  P.VC  and  then  continued  with  other  materials. 
The  fact  that  four  separate  computer"  :prcgrams  are  available  to  solve  the  disk 
flow  problem  provides  a  wide  range  of  many  materials  and  mod  if  tea  ti  oris. .  For 
rigid  PVC  degradation  kinetics  mechanism  should  definitely  be  incorporated  1  rt 
the  simulation.  In  thermosetting  resins  too,  the  increase  in  viscosity  due  to 
cross-linking  may  be  of  interest.  The  dependence  of  density,  specif ic  heat, 
thermal  conductivity  and  activation  energy,  for  flow  on  temperature  and  pressure 
will  also  have  to  be  considered. 
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APPENDIX 

Rigid  PVC  Properties  Used 

Power  law  index  h  **  .5 

Consistency  index  m  “  4.0  x  108  poise  202°C 

Pre-exponential  factor  for  the  temperature  dependence  of  viscosity 
A  =»  6.46  x  1(T8,  -  - 

Activation  energy  for  flow  AE  «  27,8  keal/g  mole 
Density  (melt)  0  »  I. 3  gm/cc 


Heat  capacity  Cp  *  ti.k 3  eal/g°K 

Thermal  conductivity  k  -  2.3  x  10  **  cal/g-secV, 

Temperature  below  which  flow  can  be  ignored  T.  *  150°C  (net  a  material 
property)  h 
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FIG.  I  THE  DISK" SHAPED  MOLD  CAVITY 
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